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Abstract 

The transition form factor F^ n (Q 2 ), decay width r(-7r° — ► 77) and charge 
form factor of pion are calculated within the model of induced nonlocal 
quark currents based on the assumption that the nonperturbative QCD 
vacuum can be characterized by a homogeneous (anti-)self-dual gluon field. 
It is shown that the interaction of the quark spin with the vacuum gluon 
field, being responsible for the chiral symmetry breaking and the spectrum 
of light mesons, can also play the decisive role in forming the form factor 
F 1 - K {Q 2 ) and decay width r(7r° — > 77). An asymptotic behavior of quark 
loops in the presence of the background gluon field for large Q 2 is discussed. 



I. INTRODUCTION 



Pion electromagnetic properties were intensively studied within different theoretical 
approaches, such as perturbative QCD [|]], the factorization theorems and QCD sum 
rules @- 0, approaches based on the Bethe-Salpeter equation fXl]- []13] , the instanton liq- 



uid model [14], relativistic quark models ||U, nonlocal quark model with confinement [1(] 



and some others. The reason for this permanent interest is quite understandable: a variety 
of general features of strong interactions clearly manifest themselves in pion physics. First 
of all this relates to the mechanisms of chiral symmetry breaking in QCD and general 
behaviour of electromagnetic form factors of hadrons. The charged form factor played the 
particularly important role for explanation of the quark counting rules |0| and factoriza- 
tion hypothesis within QCD (for review see, e.g. [0). The two-photon decay of neutral 
pion and the transition form factor relates to the triangle anomaly and is important for 
testing the models of QCD vacuum [0,§J. 

The present paper is devoted to description of electromagnetic properties of pions 
within the model of induced nonlocal quark currents developed in our previous pa- 
pers [T^JT^l • Namely, we calculate the charge and transition form factors and two-photon 
decay constant of pion. The model is based on the assumption that the (anti-) self-dual 



homogeneous field [[18H 21 



B^x) = nB^Xv, n = A 3 cos£ + A 8 sin£, 

= ±B^ U , B^ p B pfl = —B 2 S^ V (1) 

can be considered as a dominating gluon configuration in the QCD vacuum. In other 
words, it is assumed that the effective potential for QCD has a minimum at nonzero 
strength of the background field |19| , p2| . 



Vacuum field ([I]) leads to spontaneous violation of a range of symmetries such as CP, 
colour and 0(3). A satisfactory restoration of these symmetries at the hadronic scale 
assumes an existence of domain structures in the vacuum. In a given domain the vacuum 
field has a specific direction and is either self-dual or anti-self-dual, but this is uncorrelated 
with the specific realization of Eq.(|T]) in another domain. The idea of domains in the 
QCD vacuum was discussed in application to various homogeneous fields [|18]|23|-|25| . In 



complete theory the domain walls should be describable by some appropriate classical 
solitonic solutions of equations of motion. In the effective model under consideration this 
idea is realized as a prescription that different quark loops (namely, those separated by the 
meson lines) in a diagram must be averaged over different configurations of the vacuum 
field ([!]) independently of each other. 
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A consideration of both the quark-gluon dynamics in the external field and its 
manifestations in meson properties indicates that this field provides for a picture of con- 
finement and chiral symmetry breaking which is unexpected but looks quite self-consistent 
and agrees with meson phenomenology. 

Being taken into account nonperturbatively, the vacuum field under consideration 
changes analytical properties of the quark and gluon propagators drastically and makes 
QCD to be a nonlocal quantum field theory. The propagators of color charged fields 
in the presence of field (p]) are entire analytical functions in the complex momentum 



plane [jTU[ IE ,p2 which means an absence of quarks and gluons in the observable spec- 
trum of hadrons. The propagators are modified essentially in the infrared region and show 
the usual behaviour at short distances. Consideration of the Wilson loop in the presence 
of this field shows that the Wilson criterion is satisfied [^6j with an oscillator potential 
between heavy quarks. The potential arises effectively due to an interaction of the charges 
with vacuum field but not by virtue of gluon exchange. 

Due to the (anti-) self-duality and homogeneity of the field ([!]) the operator j^d^—iB^) 
has an infinite set of zero eigenmodes whose contribution to the quark propagator leads 
to the relation between densities of quark and gluon condensates: 

£ J™ n m /(9/(a0?/(*)) = _JV Fo3' ( 2 ) 

/=i 

which is valid to all loop orders in QCD. Equation (||D indicates a non-Goldstone mecha- 
nism of chiral SU^Np) x SU^(Np) symmetry breaking (for details see [f22fl ). This way of 
violation of chiral symmetry is unusual but, as it will be discussed further below, is not 
in conflict with the observable meson properties |[17|| . 

Manifestations of this quark-gluon dynamics in the spectrum of collective mesonic 
excitations were studied in [|T6|JT7[] within the model of induced nonlocal quark currents. 
Technically, the model is based on the bosonization of one-gluon exchange interaction of 
quark currents in the presence of the homogeneous (anti-)self-dual vacuum gluon field ([!]). 

The vacuum field determines basic properties of meson spectrum in the way that is 
completely consistent with experimental data. Namely, due to the zero modes in the 
quark propagator the masses of light pseudoscalar and vector mesons are strongly split, 
scalar and axial mesons as a simple qq systems in the ground state are absent in the 
spectrum (there are no scalar and axial analogies of pions and p-mesons). Zero modes 
affect drastically the pion weak decay and determine the correct value of constant /„.. An 
entireness of quark and gluon propagators results in the Regge character of the spectrum 
of orbital and radial excitations of light mesons. Moreover, scalar and axial mesons 
with appropriate masses (e.g., a\ meson) appear in the hyper-fine structure of orbital 
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excitations of vector mesons. The mass of heavy quarkonium tends to be equal to sum of 
the masses of quarks, the heavy-light meson mass approaches a mass of a heavy quark, and 
the weak decay constant for pseudoscalar heavy-light mesons has asymptotic behaviour 
1 / y/rriQ, which is consistent with the Isgur-Wise symmetry. Quantitatively, the masses 
and decay constants of mesons from all different regions of the spectrum are described 
within ten percent inaccuracy. These different phenomena are displayed with the minimal 
set of parameters: gauge coupling constant, strength of the vacuum field and the quark 
masses. 

This paper proceeds a systematic analysis of possible manifestations of the vacuum 
field ([[]) in meson properties. We intend to give a kind of unified description of a wide 
class of static (masses, decay widths) and dynamical (form factors) characteristics of 
different mesons with the set of parameters, that is minimal for QCD. Using the values 
of parameters fitted from the meson spectrum we got the electromagnetic form factors of 
pions to within 15% inaccuracy. . 

Our main goal is to study the role of quark zero modes, induced by an interaction of 
quark spin with the vacuum field, in the two-photon decay and transition form factor of 
7T°. The main result is an observation that the vacuum field ([!]) can be responsible for 
both of them. Via zero modes the vacuum field affects in a crucial manner the form factor 
and decay constant. In particular, we show that the spin- field interaction generates the 
triangle anomaly in the pion decay to two photons and leads to quantitatively satisfactory 
description of the transition form factor and decay width. 

The main effect of the vacuum field under consideration in the charge form factor is an 
increasing of the contribution of triangle diagram to the form factor at moderately large 
transfer momentum Q 2 . This effect is due to the presence of the background field both in 
the quark propagators and nonlocal meson-quark vertices, which causes a specific interplay 
of translation and color gauge invariance in the quark loops. Namely, the translation 
x — > x + a causes a shift in the field B^ v x v — > B^ y x v + B^ u a V) which can be compensated 
by an appropriate gauge transformation. Therefore, only gauge invariant quantities turn 
out to be translation invariant. This results in nonconservation of energy- momentum 
in the separate vertices of Feynman diagrams, but the conservation is still held for a 
whole diagram if it is gauge invariant. This peculiarity changes an asymptotic behaviour 
of some Feynman diagrams at large momentum transfer. This is particularly relevant 
to the triangle diagram for the charge form factor. Usually the asymptotic behaviour 
~ (1/Q 2 )^ 2 of this diagram is a direct consequence of the behaviour of a meson-quark 



vertex and quark propagator at short distances |fLl| . The presence of the vacuum field 
changes the asymptotics cardinally, and it becomes ~ l/(Q 2 ) 1+m V B , where m q is the 
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quark mass and B is the strength of the vacuum field. This improves the form factor at 
intermediate and moderately large Q 2 . However, the absolute asymptotics of the form 
factor within the model under consideration comes from the so-called bubble diagrams. 
These diagrams can be treated as describing the hard rescattering of quarks inside a pion 
via virtual gluon exchange. This mechanism is in agreement with the general analysis 



within quark counting rules and perturbative QCD 0,pl| . 

In the next section we review shortly the model and introduce the electromagnetic 
interactions into its structure. Section III is devoted to calculation of transition form 
factor, where we analyze the role of quark zero modes. Charged form factor of pion is 
calculated and the manifestations of the vacuum field in asymptotics of Feynman diagrams 
at large momentum transfer is discussed in sect. IV. Several Appendices contain details 
of calculations. 

II. THE MODEL OF INDUCED NONLOCAL QUARK CURRENTS 
A. Basic Approximations and Notation 

The functional integral for Euclidean QCD in the presence of the vacuum field (|]) can 



be written in the form 22 



N 



Z = N ( do B f Dn G (G,B) f \{Dq f Dq f 

D/2 G (G, B) = DGA FP [G, B}5 [V(B)G]exp y d 4 xC YM [G + B] j , (3) 
G = lfl G„, V = 7/i V M , V M = <9 M -^. 

The vacuum field B^ (gB = B) is characterized by angle £ and two spherical angles (ip, 9) 
defining the directions in the color and Euclidean spaces respectively. The measure of 
integration das is defined as 

± V ' 

where the symbol J2± denotes averaging over the self- and anti-self-dual configurations. 

The functional space Q contains the quantum gauge fields G^ vanishing at the space- 
time infinity. The fermionic functional integral spans the Grassmann algebra T of square 
integrable quark fields. An anti-hermitean representation for Euclidean 7-matrices is 
adopted. 



5 



The functional integral (|3|) contains the intrinsic dimensionful quantity B, which is a 
gauge and renormalization invariant strength of the vacuum field, which by assumption 



minimizes the QCD effective potential (e.g., see ]19 , 18| , 20] ) . This field strength (or gluon 
condensate) provides the natural reference scale for running quark masses and 
gauge coupling constant a s (fi). Therefore, the strength of the vacuum field B, the quark 
masses and coupling constant at the scale /i = \/B can be considered as the physical 
(intrinsic) parameters of QCD in representation (|3|). The strength B can be related to 
the fundamental scale Aqcd |T9|J2^] . Values of the parameters have to be extracted from 
the analysis of hadron spectrum. To be able to do this one needs to rewrite representation 
(|3|) in terms of composite hadron fields, describing collective colorless excitations in the 
QCD vacuum. For meson fields this program has been realized in papers pB|JT7| , where 
the model of bosonization, referred below as the model of induced nonlocal quark currents, 
has been developed, and the masses and weak decay constants of mesons from the different 
regions of meson spectrum have been calculated. The values of parameters and the results 
for meson masses are summarized in Table | and Fig. 1. 

A derivation of the model goes through the following steps. Integrating over gluon 
fields G one can rewrite Eq. ([3]) in the form 

N f 

Z = N [ da B [ \{Dq f Dq f 

exp | J d 4 xq f (x) (iV - ra/) q f (x) + ^ £ n j , (5) 

Ln = nil d " Vl -S d4y ^ y ^ ■ ■■■jZ(yn)G%:X(yi,..,y n I b), 

Np 

j;(y) = J2qf( x h, ta( if( x )- ( 6 ) 

/ 

The functions Gf}' "f™ are the exact n-point gluon Green's functions in the external field 
Bl- 

The model of induced nonlocal quark currents is based on the assumption that the four- 
quark interaction L 2 in Eq. fl5|) plays the main role in the formation of mesonic qq collective 
modes, but the rest of terms L n (n > 2) can be omitted in the first approximation. Thus, 
we consider Eq. ([|) with the quark-quark interaction truncated up to the term L 2 

N f 

T n„ I f j4„- /„\ f.-vr ~_ "\ ~ f~\ i 7 

^2 



Z = N J das J Y\_DqfDq~f exp I^J d 4 xqj(x) (iV — m^j qf{x) + L 2 

If d'xdS j;(x)G a Jl(x,y | B)jl(y). (7) 



9 2 
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Gluon two-point function G^(x,y \ B) is approximated by the gluon propagator in the 
external field fll]), so that the vacuum field is taken into account exactly but the radiative 
corrections are neglected in the model. The subsequent steps are straightforward and 
consist in a standard Fiertz transformation, decomposition of L2 into an infinite series of 
current- current interaction terms and bosonization of the effective four-quark interactions. 



For details we refer to papers [|16| , |17jl . 

The starting point for incorporating the electromagnetic interactions into the model 
under consideration is the representation [ITj 



I+SJ 

E E 

aJin j=\£-sj\ 



2A2 M.Un 



d x 



■j-ajinj ! 



X 



(8) 



where 



a = 0,l,...,N f , J = S,P,V,A, 
sp = s s = 0, s v = s A = 1, 



0,1,..., 71 = 0, 1,..., 



G 



Jin 



Cs = C P = 1/9, C v = C A = l/U 



(9) 



2 e n\(£ + n)\ 

The currents X a,Un i are nonlocal, carry a complete set of quantum numbers (isospin a, 
spin-parity in the ground state J, orbital £ and total angular j momenta, and a radial 
number n) and has the following form 

^ut = J"%e fOT J = S,P, £>0 (j=£), 
X^om = juon for J = ViAi i = ( J = l), 



1 



,3=* 



O.H1--HI 



(£+1)2 ' a/Jti—Hi 
1 y4 [ jbJin _ jbJln 



7 
>->< 



bJin 



TbJln 



j=£+l. 



(10) 



- -M 7 

Symbol V a Hi...Ht denotes a cyclic permutation of the indices (a/xi...//^). The currents J 
in Eq. ( [I0"D read 



6J£n 



,H\---HV 



7° 

irbj £n 
' 1 



q f [x) 



ybJln 



a, m. ..^e 



X 



ff 



F ni {4s) 



M b ff ,T J 

ff OL 



dtt l+n e s \ 




A 2 



7^) 
in—to 



1 V//' (g 
A 




V//'= 6/ V -£// V, X7=d-iB, \7=d+iB, = m f /{m f + m f ), 



IF 



«75, 



y — 7 



7s7 
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The polynomials a,{ x ) are the irreducible tensors of 0(4). The doubled brackets in 
Eq. ( PD mean that the covariant derivatives commute inside these brackets. For the sake 
of simplicity the scale A 2 = y/3B is introduced and a particular direction n a = 5 a8 of the 
vacuum field in the color space is fixed, so that: 



v = diag(l/3,l/3,2/3). 



T21 



Further details concerning above representations can be found in papers PqJTi 

Form-factors F n e(s) appeared as a result of decomposition of bilocal quark currents 
over the orthogonal complete set of polynomials for which gluon propagator plays the role 
of a weight function. The are entire analytical functions in the complex s-plane, which is 
a manifestation of the gluon confinement. Now an interaction between quarks is expressed 
in terms of the nonlocal quark currents, which are elementary currents of the system in 
the sense of classification over quantum numbers. 



B. Electromagnetic Interactions 



Interaction of quarks with electromagnetic field can be introduced into represen- 
tation (0) with L 2 written in the form of Eq. @ by means of the minimal substitution 
(see also [§]) 



V^V —iefA(x) 
V^V +ie f A(x) 



(13) 



both in the free quark Lagrangian and vertices V a ji n j- Here e/ = eQf, and Q is the charge 
matrix. 

Equation (0) takes the form (TV = {aJinj}) 



(14) 



Z[A] = j da B J DqDqexp^J d xq f (x)[i^ V M -rrif + e/7 M A j[i (a;)]g/( 
+ E^ G ^ / ^x\l M (x I A)-TrV M {x | 0)S(x,x | 0)| 2 }, 



where the currents Xj^{x \ A) include electromagnetic field A through the vertex functions 



V bJtn (x | A) = M b T J U F 7 



b r J 



/«2 

D (■'•) \ T(£) I ID Cr 



D ff ,= £/ 



V +ie f A^(x) 



\ 



A 2 



% A 




V -ie f A^(x) 



(15) 
(16) 



One can check that under the gauge transformation 
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q f (x)^q' f (x)=e ie f a ^q f (x), 
q f (x) -> q' f (x) = q f (x)e-^\ 

A^x) -> = A^z) + <9 M a(x) (17) 

the currents Tjj are transformed with the electric charge (e/ — ey) corresponding to a 
meson composed of qj and g/': 

X N {x | A) -> JV(a; | A) = ^~ e f'^Z M {x \ A), 

l^{x | A) -> 2j(a; | A) = J^(x | A)e- <(e '- e /' )a(a:) . ( 18 ) 

The current-current interaction lj^(x)Xj^(x) is invariant under U(l) transformations (|17|) . 
By means of the standard bosonization procedure |l6| , p!7| applied to Eq. (|14]) the 



functional integral Z[A] can be represented in terms of the composite meson fields $/v" 
Z[A] =N J n^Arexp [±JJ d 4 xd 4 y^ M (x) [(□ - M 2 -) 6{x - y) 
- h 2 M U^x - y)} $„(y) + I int [$ | A]} , (19) 



/int = - / d*xh N <b N {x) [T N {x | A) - rV(a; | 0)] 

-- j d^xi J d 4 x 2 h A fh A f l § A f(x 1 ) [T A f A f / (x 1 ,x 2 \ A) - S A f A f / U Af (x 1 - x 2 )} <&u>{x 2 ) 
lr r m 

~ d 4 X!... / d 4 x m Y[ hM k $M k (x k )rM 1 .„jv m {x 1 , ...,x m | A), 

m=3 m J J k=l 

^ATi..jf m = J da B Tr {V^ixx | A)S(x 1 ,x 2 \ A)...V Afrn (x m | A)S(x m ,x 1 \ A)} . (20) 

At zero electromagnetic field the effective meson action in Eq. flT9"|) coincides with the 
action derived in |T7]]. Meson masses are calculated by solving the equations 

A 2 + G\U M {-Ml) = 0, 

where Ilj^(p 2 ) is a diagonal part of the Fourier transform of the two-point function IV^v 
at zero electromagnetic field. 

The meson-quark coupling constants are defined by the relations: 

^ = i/nV(p 2 )| p2= _M^. (21) 

The integral over directions of the field and the sum over self- and anti-self-dual con- 
figurations has been transferred to the exponent in Eq. ( |T§D and included into the vertices 
T, Eq. (p0[) . This assumes that the vacuum field configurations in the different (sepa- 
rated by meson lines) quark loops are considered as independent from each other. Such 
a prescription should reflect effectively the domain structures in the vacuum. 
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To get various meson-photon amplitudes one has to decompose the vertex functionals 
~^Ni...N m into a series over the electromagnetic field A, as is shown schematically in Fig. 2. 
This is achieved by a decomposition of the quark propagators, 

Sf(x,y\A) = S f (x,y) 

+ ^e n / d A z l ... J d 4 z n Sf(x,z 1 )-f fll A fll (z 1 )...-f fln A fMn (z n )Sf(z n ,y), (22) 

n=l 

and vertices 

V M {x\A) = V M {x) 

+ £ e n J d A Zl ... j c^„41...m>> *i, -. *04* (zi)...^^). (23) 



n=l 



A regular procedure for calculation of the vertices Vfjr/^ „ is described in Appendix A. 

The lowest vertex for charged pion, the case of particular importance for further 
calculations, comes from the current J~ w +, Eqs. (|TTD and (]I5D , 



X+ix) = d{x)i^F 00 




u(x) = d{x)i^^ 




U\X) 



+i7 5 e J d 4 "yd(x)A fl (y)T fl (x,y)u(x) + 0(e 2 ), 



(24) 



where 



IV(z,y) 



i i 



(2tt)< 



o o 



4A 2 



2z V M (x) - <^ 



x exp 



14A 2 



V (x) + (3(2i V (x) - q)q 



(25) 



The first term in Eq. fl24j) is the nonlocal quark current in the absence of field but 
the second term describes an interaction of photon with a quark inside pion which is 
characterized by the form factor T^(x,y) f8|. 

The quark propagator in the homogeneous (anti-) self-dual gluon field Sf(x,y) is a 
solution to equation 



(^V M - m f )S f (x, y) = -6(x-y), 
S f (x,y) = e^ XliB ^ Vv Hf(x - y), 



Hf(z) = (ij^Vn + mf) f-V + m f - a af3 B a(3 ) 5( 



-i 



(26) 



The term <J a pB a p in Eq. (|26[) is particularly important. It describes an interaction of a 
quark spin with the vacuum field and is responsible for the quark zero modes. Fourier 
transform of the translation invariant part Hf reads [17 , 22 
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1 / 4r 
** f ^ = 2vA^ I dse ~*^ S (1+7) ^ A ][Vala ± is TBTa/«j9P/9 



/ 1 + i 

+m f P± + P T - - -7a/a/37/3 



1-s 2 2 '" J " pip l-s 2 

A 8 

vA 2 



(27) 



where P± — ( I ± 75) /2 and the upper (lower) sign corresponds to (anti-)self-dual config- 
uration. 

Propagator fl27D is an entire analytical function in the complex momentum plane, 
which is treated as quark confinement. The term a a pI3 a ^ in Eq. (p6|) is particularly 
important. It describes an interaction of a quark spin with the vacuum field and is 
responsible for the quark zero modes. Contribution of zero modes to the propagator is 
seen in Eq. (|27|) as a singilarity of the integrand at s — 1 which is integrable unless rrif = 0. 

Below we show that the spin-field interaction is of crucial importance for the transition 
form factor of neutral pion. The presence of the vacuum field in the phase factor in quark 
propagator and in covariant derivatives in the vertices turns out to be important for the 
charge form factor. 



III. DECAY 7T° -> 77 AND 7*vr u -» 7 TRANSITION FORM FACTOR 

A vertex relevant to the interaction of neutral pion and two photons with momenta p, 
ki and k 2 , correspondingly, has the following general structure 

T^ip, K h) = iS^ih + k 2 - py^k^Mp 2 , kf, k 2 2 ), (28) 

where T^^p 2 , k\, k\) is a scalar function. It should be remembered that we have started 
with the Euclidean formulation of QCD, and fcj, k% and p are Euclidean momenta. 

For 7*7T° — > 7 transition the final photon 7 is on the mass shell k\ = 0, whereas 
k\ = Q 2 > for the virtual photon 7*. The transition form factor F ln (Q 2 ) is then defined 

as 

F^(Q 2 )=T 7r (-M 2 ,Q 2 ,0). (29) 
The two-photon decay coupling constant 

g^ = T w (-M 2 , 0,0) (30) 

and the decay width 
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IV - 77) = 4« 2 M^ 77 . (31) 

correspond to kinematics with on-shell photons and pion: k\ = k\ = and p 2 = —M 2 . 

As follows from the effective action in Eq. (|20|) , the one loop approximation of the 
vertex (p8| ) is described by the triangle diagram shown in Fig. 3. Other one loop diagrams 
with photon lines attached to the meson-quark vertex are identically equal to zero due to 
the parity conservation. 

As has been mentioned above, the interaction of quark spin with the homogeneous 
vacuum field generates an infinite number of zero modes of Dirac operator and leads to 
a nonzero quark condensate density, which indicates a breakdown of chiral symmetry by 
the vacuum gluon field . A contribution of zero modes to the polarization diagrams of 
light mesons is responsible for mass splitting between vector and pseudoscalar mesons and 
smallness of pion mass [17]. The main goal of this section is to show that the zero modes 



generate the triangle anomaly and play the decisive role in the two-photon decay of 7r° 
and transition form factor F 17T . To demonstrate this qualitatively it suffices to consider 
the triangle diagram (Fig. 3) in the limit of vanishing quark masses rrid = m u — > 0. 

The quark propagator in the external (anti-)self-dual field (0) has the following stan- 
dard representation in terms of the matrix elements of the projection operators V n onto 
the subspaces corresponding to different eigen numbers A n of Dirac operator in the exter- 
nal field HI! 

One can separate the contribution of the zero eigenmodes and normal modes to the 
propagator 

S(x,y) = S'(x,y) + S (x,y), 

S'{x, y) =ii/ x A{x, y)P ± + A{x, y) iV y P T + 0(m), (32) 
\7=d-iB, \7=d+iB, 
S (x,y) = V (x,y)/m. 

Here A(x,y) is the scalar massless propagator in the background field (|1|), and P± = 
(1 ± 75)/2. The projector onto the zero mode subspace looks as 

n 2 B 2 

V = — r f(x,y)P T ^ (33) 

7T Z 

S± = - (1 ± T,jbj) , Sj = -e ijk cr jk , a jk = [7^, -y k ]/2i, 
bj — Bj/B, Bi = —-Eij k Bj k , i,j,k = 1,2,3 
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f(x, y) = exp l--Vr)?B(x - y) 2 + inx^B^yJi , (34) 
d 4 zV (x, z)V (z, y) = V (x, y). 



The projector Vo(x, y) contains the projection matrices P_ (self-dual field) or P + (anti-self- 
dual field) and S_ {n > 0) or E + (n < 0). The matrix S + (£_) can be seen as the projector 
onto the quark state with the spin orientated along (against) the chromomagnetic field 
By Now consider the n-point quark loop 

TriySoOri, x 2 ) + 5"(xi, x 2 ))r 2 (S'o(a; 2 , x 3 ) + S'(x 2 , x 3 )) . . . T n (S (x n , x^ + S'(x n , a?i)), 

(35) 



where are some Dirac matrices. Chiral structure of zero mode part Vq, Eq. (|33|), and 
normal mode part S', Eq. (]32f) , suggests that the loops with all vector vertices are regular 



in the massless limit (for details see f22[), while diagrams with one pseudoscalar and n — 1 
vector vertices are singular and behaves as l/m u . The second kind is just the case of 
Fig. 3. 

Thus the triangle diagram for the amplitude ( p8|) is proportional to l/m u as m u — > 0. 
However, the amplitude contains also the pion-quark coupling constant h n defined by the 
pion polarization function according to Eq. (pT|). The polarization function is simply a 
two-point quark loop of the form fl35|) with pseudoscalar vertices. In the massless limit 
it diverges as l/m 2 and, hence, coupling constant h n oc m u . Thus the massless limit of 
the form factor, that is a product of coupling constant h n and the quark loop, is nonzero: 
lim^^o F-y-w 7^ 0. This anomaly is determined completely by a contribution of quark zero 
modes to triangle diagram and coupling constant. 

Let us consider this important peculiarity in more details. Using the proper time 
representation of the quark propagator (27) and meson-quark vertex (|TT|) , evaluating the 
trace of Dirac matrices, calculating the loop momentum integrals and averaging over 
different configurations of the vacuum field, we represent the form factor F 11T (Q 2 ) as an 
integral over proper times t and si, s 2 , S3, corresponding to the vertex and propagators, 
respectively (see Fig. 3): 



1 1 



A2V2?r 

1 ^ 



si\ fl - s 2 \ (I - S3 



1 + S1/ vi + s 2 y V1 + S3 



ml/4v 



Tfl r 

x E T ^ 2 <f l (M2,Q 2 ;s,t)exp[M^(s,t)-QV(^^l, (36) 

i=l,2,3 

0(s, t) = (2sis 3 + vt(si + s 3 - s 2 (l + SiS 3 )))/4t>x, 
<p(s, t) = s 2 (s 3 + vt(l + s 1 s 3 ))/2vx, 
X = 2vt(l + sis 2 + S1S3 + s 2 s 3 ) + (si + s 2 + s 3 + sis 2 s 3 ). 
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Here and below we use the shorthand notation for dimensionless ratios: Q 2 = Q 2 /A 2 , 
m u = m u /A, M n = M n /A. The symbol Tr„ means summation over the elements of 
the diagonal matrix v written in Eq. (|12|) . The details of calculation of F 17T (Q 2 ) and an 
explicit form of functions &i(M 2 , Q 2 \ s, t) can be found in Appendix B. 

The singularities (1 — of the integrand in Eq. fl36|) at Sj — > 1 appear from the 
zero mode contribution to the quark propagator (see the second line in Eq. (^)). These 
singularities lead to the 1 /m^-dependence of the integral in Eq. (|36"D in the limit m u <C 1: 



f,AQ 2 ) = —HQ 2 )- (37) 

Here / does not depend on m u . In the massless limit pion polarization function 11^ (— M 2 ) 



looks as 17 



mi n l Mi 




Ml 



exp < > — exp 



3v(l + 2v) 

hence the effective coupling constant h n in Eq. (|37D behaves as 



(3* 



K = 1/ Vn;(-M2) ~ m u , (39) 
and F 17T (Q 2 ) does not vanish as m u — > but approaches a constant value. The same is 



valid for decay width (|31"|). 

Numerical results for the form factor and decay width are represented in Fig. 4 and 
Table II. The model parameters are given in Table I. It should be remembered that we 



simply use the values of parameters fixed from the description of the meson spectrum [|T7 



The solid curve in Fig. 4 corresponds to Eq. (pq). The radius for 7*71- — > 7 transition 
defined as 

<I -" > V 7I (0)' 

is equal to .57 fm, that have to be compared with r^ p = .65 ± .03 [pCfl . 

Just to illustrate the crucial role of the quark zero modes we show the form factor 
calculated with zero mode contribution eliminated from the quark propagator (the long 
dashed line). Zero modes drastically affect also the two-photon decay constant g irn and 
decay width r(7r° — > 77), as is seen from Table II. 

Above consideration supports the picture of chiral symmetry breaking due to the 
fermion zero modes induced by the homogeneous (anti-)self-dual vacuum gluon field, 
which has been developed in our previous papers fl7| and |22| . We can conclude that 



within the model of induced nonlocal currents the experimental data for 7*71° — > 7 tran- 
sition form factor and two photon decay constant can be explained by the same reason as 
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a smallness of pion mass, splitting of the masses of vector and pseudoscalar light mesons 
and weak decay constants of pions and kaons. A general physical reason is an interaction 
of a quark spin with the vacuum homogeneous gluon field. This spin-field interaction 
appears to be a dominating effect in the above discussed phenomena. 

It is appropriate to mention here, that, possibly, there is another, completely inde- 
pendent of our considerations, manifestation of an interaction of the quark spin with a 
long-range gluon field in the QCD vacuum. As is reported in paper |HJ the experimen- 
tally observed sign of the jet handedness correlation can be naturally understood if the jet 
fragmentation occurs in the background of vacuum gluon field that is (almost) homoge- 
neous within some characteristic region. The spin-field interaction plays here the leading, 
qualitatively important, role. 

The asymptotic behaviour of the Feynman diagrams in the limit of large momentum 
transfer is an additional point where the homogeneous vacuum field could be seen as a 
relevant effect. We will pay more attention to this in the next section, where the charge 
form factor of pion is considered. However, just for comparison, it is advantageous to 
calculate the asymptotic form of the transition form factor. 

A behaviour of the triangle diagram for the transition form factor in the limit Q 2 ^> A 2 
can be easily estimated. Equation (|36f) can be rewritten in the form 

F^(Q 2 /A 2 ) = -J ... J dtds 1 ds 2 ds 3 ^(Q 2 s 2 /A 2 ,M 7r /A,m u /A;t,s u s 3 ) 



" Q 2 

xexp -—(p(s,t) 
<f(si, s 2 , s 3 , t) = s 2 (s 3 + vt(l + Sis 3 ))/2vx{si, s 2 , s 3 , t), 

which just underlines that the integrand $ depends on Q 2 and s 2 in the combination Q 2 s 2 , 
where the variable s 2 corresponds to the quark propagator situated in Fig. 3 between two 
electromagnetic vertices. Here we have restored the dimensional notation for the masses 
and momentum Q. Function \ is given in Eq. fl3"6|). For any fixed values of t, si, S3, 
the function ip is increasing in s 2 and gets the lowest value at the point s 2 = 0. This 
corresponds to the ultraviolet regime in the space of integration variables, that is usual 
for large momentum asymptotics 0. An integral over s 2 can be evaluated by the Laplace 
method and, for Q 2 ^> A 2 , we arrive at relation 

Q 2 F^(Q 2 /A 2 ) = CA + 0(A 2 /Q 2 ) « 0.2Gev + 0(A 2 /Q 2 ), 
11 00 

C = Tr v f ... J dtds x ds 3 J ds 2 $(s 2 , M w /A, m u /A; t, Si, s 3 ) 

00 
x exp [-S2O3 + vt(l + s 1 s 3 ))/2vx(s 1 , 0, s 3 , £)] . 
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This result has to be compared with the Brodsky-Lepage limit [1] 

Q 2 F 17T (Q 2 ) -> 2F n = .186Gev. 

We see that the asymptotic regime in the form factor is realized in the usual way consistent 
with the analysis within factorization hypothesis and QCD sum rule approaches ||. 



IV. THE PION CHARGE FORM FACTOR 



According to the effective action (|20D , the one-loop amplitude for the processes 
7r ± 7* — > 7r =l= is described by the triangle and bubble diagrams shown in Figs. 5a and 
5b correspondingly. It has the following structure 



A"(fc ls k 2 , q) = $< 4 > (h -k 2 + q) \{k x + k 2 fF 1 (k 2 1 , k 2 2 , q 2 ) + q,F 2 (kl k 2 2 , q 2 



(40) 



where k\ and k 2 are the pion momenta and q is a momentum of virtual photon 7*. 
As is known, within the minimal substitution scheme of introduction of electromagnetic 
interactions the vertex A M given by a sum of triangle and bubble diagrams satisfies the 
Ward-Takahashi identity [|8]|§. 

The pion charge form factor F n (Q 2 ) is defined by the relation 

F n (Q 2 ) = F n A (Q 2 ) + F:(Q 2 ), 



F^-Ml-MlQ 2 ), F: = F°(-M 2 ,-M 2 ,Q 2 ). 



The details of calculation of contributions of triangle F^{Q 2 ) and bubble F°{Q 2 ) diagrams 
to the form factor are relegated to Appendix C. The result expressed in the form of the 
proper time integrals reads: 



—2-Tr 

4tt 2 " 



dtidt 2 ds\ds 2 ds 3 



si 



S 2 



S 3 



1 + siJ VI + s 2 J VI + s 3 , 
Q 2 ^{Ml Q 2 - s, t) + m 2 u 1> 2 (M 2 , Q 2 - s, t) + M 2 <S> 3 (M 2 , Q 2 - s, t) 







+^(M 2 ,Q 2 ;s,t) exp M 2 <f>{s,t) - Q 2 <p{s,t) 



(41) 



(j) = (2sis 2 + 2sis 3 + v(tt + t 2 )(s 1 + s 2 + s 3 + sis 2 s 3 ))/4t;x, 

X = 2v(t 1 + i 2 )(l + sis 2 + S1S3 + s 2 s 3 ) + (1 + Av\t 2 )(s 1 + s 2 + s 3 + S1S2S3), 

f = [s 2 s 3 + vtis 2 (l + S1S3) + vt 2 s 3 (l + sisz) + vHxt^l + sis 2 + sis 3 + s 2 s 3 )]/2vx, 

for triangle diagram, and 
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hi 



— Tr 



4?r 2 

2 



dtidt2(is\ds-2,d(3 



1 + sJ U + sJ U + S3 



1 m£/4i> 



-$°(M 2 , Q 2 - 0, s, t) + <&°(M 2 , Q 2 - (3, s, t) 



x exp 

0° 

X° 



M 2 <f(s,t)-QV(/W) 

[2s 1 s 2 + v(t 1 + t 2 )(s 1 + s 2 )]/Av X 
vt 2 /3 



2v X ° 



[s! + vti(l + sis 2 ) + ut 2 (l -/?)[! + sis 2 + 2vt 1 (s 1 + s 2 ) 



2u(*i + t 2 )(l + s x s 2 ) + (1 + 4w 2 M 2 )(si + s 2 ), (42) 

for bubble diagram. We have used the shorthand dimensionless notation for momentum 
Q and masses. Functions $j and $° are written in Appendix C. 

For the parameter values given in Table 1 the charge form factor defined by Eqs. (41) 
and fl4"2"D is plotted in Fig. 6 by the solid line. The electromagnetic radius takes the value: 

•^(0) _ 



(<) = -6 



F^oy 



.524 fm. 



One sees that agreement with experimental data for the form factor and radius r° xp = 
.656 fm is quite satisfactory. 

An improvement of the radius and form factor at small Q 2 can come from the diagram 
with intermediate p-meson (Fig. 7). It is known @JT^| that its contribution to F n (Q 2 ) 
can be important in the region Q 2 < 5Gev 2 . An estimation of the diagram with p-meson 
within the model under consideration also indicates a diminishing of the form factor, 
which is maximal (about 6%) in the region of Q 2 ~ 2Gev 2 . 

Numerically the contribution of the bubble diagrams is very small. It is maximal at 
Q 2 ~ 2Gev 2 and is of order 10~ 3 . Thus within the model under consideration the triangle 
diagram gives the main contribution to the form factor for the values of Q 2 shown in 
Fig. 6. One sees that the calculated form factor (the solid line) smoothly approaches the 
experimental fit (dashed line) at large Q 2 . This behaviour seems unexpected. As is known 
from the studies, based on the effective meson-quark quantum field models [fsHTT|, the 
triangle diagram should decay stronger than 1/Q 2 - the asymptotics of the experimental 
fit. In our case, a naive estimation based on the ultraviolet behaviour of the quark 
and vertex (|Tl~D gives (Q 2 )~ 2 . However, the homogeneous vacuum gluon 



propagator 

field changes the asymptotics of the triangle diagram cardinally. 

To demonstrate this let us consider the function in the limit Q 2 ^> A 2 . It is 

convenient to rewrite Eq. (f4T|) as 

1 1 



hi 
^Tr, 



4vr 2 



dt\dt 2 ds\ds 2 dsz 



S2 



l + s 2 



1 mi/Av 
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x$ A (Q 2 , M 2 , ml; s l7 t t ) exp j-|^(g 2 , 8 A , (43) 

$ A = exp (-Q 2 ^) [g 2 $ !(M 2 , g 2 ; s, t) + ml$ 2 (Ml Q 2 ; s, t) 
+M 2 $ 3 (M 2 , g 2 ; s, t) + $ 4 (M 2 , g 2 ; s, t) 

Function Q has the form 

^ nt ,\ m?„ l + si 1 - s\ m 2 . 1 + si , t . s 
2Q Z 1 — si A\ + S\A 2 2Q Z 1 — s\ 

where A\ and A 2 are functions of t\, t 2 , s 2 , and S3 and do not depend on sj. Function 
is given in Eq. flUD , and 

tp = [sis 2 s 3 + wtis 2 (si + s 3 ) + vt 2 s 3 (s 1 + s 2 ) + vH^sx + s 2 + s 3 + sis 2 s 3 )]/2wx 

comes from the hyperbolic functions which have appeared in the integrand due to av- 
eraging of the quark loop over directions of the vacuum field. For details we refer to 
Eqs. ( |58| ) and fl6l| ) in Appendices B and C. We simply joined the exponentially increasing 
part of these hyperbolic functions with the exponent in Eq. (fHP- It should be stressed 
that the origin of such an exponentially increasing with Q 2 — > 00 terms in the integrand 



is the presence of phase factor exp^x^B^y^) in the quark propagator ( p?|) and covariant 
derivatives V = d — iB in the vertices flTTp. From the physical point of view, this simply 
means that in the presence of external field the translation invariance holds for the gauge 
invariant quantities but not for the vertices and propagators separately. As is stressed in 
Appendix C, this is the reason why the energy-momentum is not conserved in the separate 
vertices, but the conservation law is still valid for a whole diagram, if it is gauge invariant. 
One notices that the function Q has a minimum at Si = s*: 
dQ „ m 2 



> S *1 = 1-7^(^1+^2; 

2g 2 



dsx ' 1 2g 2 

d 2 n ^ 2 



> 0, 



3s 2 n=-I m\{A x + A 2 ) 2 
, , (l + s 2 )(l + s 3 ) (l + 2vt l )(l + 2vt 2 ) 

/±l + /i 2 



s 2 s 3 + vtxs 2 {\ - S3) + vt 2 s 3 (l - S 2 ) + f 2 tit 2 (l - s 2 )(l - S3)' 
Vti, h, s 2 , s 3 , g 2 s{ E [0, 1], lim si = 1, (45) 

Q 2 ^oo 

which allows to integrate over variable Si using the saddle-point approximation. The 
result is 

1 1 

Ff{Q 2 ) = Tr v J ... J dt x dt 2 ds 2 ds 3 ^ A {Q 2 , M 2 , m 2 ; s{, s 2 , s 3 , h, t 2 ) 




m 2 u \ mi/Av ^TrvmUA, + A 2 f f m 2 u / A 1 + A 2 

exp < — — - 1 — In 



2g 2 I Q 2 " I Av 
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Since functions contained in $ A (see Eq. ([y]) and Appendix C) have the following 
asymptotic form 

$ 2 exp (-Q 2 if) ~ 1/Q 2 (1 - sf) = 2/mKAi + A 2 ) = const, 

\ / ST=S, 



si=s; 

2 rf\ „, — ( /n2„ 



$1,3 exp . ~l/Q% $ 4 exp -Q 2 ip 



~ const, 

Sl=S* 



hence $ A does not depend on Q 2 in the leading order 

lim $ A (Q 2 , Ml m 2 u - s{, s 2 , s 3 , h, t 2 ) = $ A (M 2 , m 2 ; s 2 , s 3 , h, t 2 ) + 0(1/Q 2 ). 

Q — >oo 

Finally, the asymptotic formula for the triangle diagram reads 
where factor C A is independent of Q 2 



C A (M 2 /A 2 ,m 2 JA 2 ) 
(Q 2 /A 2 ) 1+ml/4vA2 



C A = hl^=^- ( ^ ] e 
*y/2^n V2A 
i i 



x 



/•rr> \ 1+mf. /2vA 2 
( Hht) e ~ml/ivA 2 

\2AJ 

J... J dt.dhds.ds^UM'/A 2 , m 2 J A 2 - s 2 , s 3 , t u t 2 ){A x + A 2 ) 1+m ^ 2 . 





Here we have returned to the dimensionful notation for the masses and momentum Q. 
For the parameter values from Table I we obtain the following result 

W"/A'> = (47) 

This asymptotic formula fits well the solid curve in Fig. 6 for Q 2 > 5Gev 2 . 

Thus the main effect of the vacuum field under consideration in the charge form factor 
is an increasing of the contribution of triangle diagram to the form factor at large Q 2 . As 
has been explained, this is clearly due to the presence of the background field both in the 
quark propagators and nonlocal meson-quark vertices, which causes the specific interplay 
of translation and color gauge invariance in the quark loops. A comparison with the large 
Q 2 behaviour of the ^77 triangle diagram, which is usual 1/Q 2 (see the previous section), 
indicates that a number of nonlocal vertices in a loop is of crucial importance. 

The contribution of bubble graph (see Fig. 5) to asymptotic behaviour is easy to derive 
if to rewrite Eq. ( P^| ) in the form 

1 1 

F°(Q 2 ) = Tr„ J ... J dt x dt 2 ds x ds 2 d$ $°(Q 2 /3, M 2 , m 2 ; s is U) exp {-Q V(/3, s l7 U)} , 
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which underlines that the preexponential factor $° depends on Q 2 in the combination 
Q 2 P, as is seen from functions $° (see Eq. ( |42"D and Appendix C). Here (3 is a proper 
time corresponding to the vertex operator Y ^ given by Eq. fl25|) . The smallest value of 
the function (p° corresponds to the point (3 — for any £i, £2, si, and s 2 : 

y?°(/3, ti, £ 2 , si, s 2 ) = /9<^ (/5, £1, £ 2 , s x , s 2 ), 



dp 



(3=0 



V9 o (0,£i,£ 2 ,Si,s 2 ) 



w£ 2 
2^ 



si + + £ 2 )(1 + sis 2 ) + 2w 2 £i£ 2 (si + 8 2 ) 



> 0. 



Therefore the leading term takes the form 



A 2 



'A 2 



F;(Q 2 /A 2 ) = C°(M 2 /A 2 ,m 2 /A 2 )— + O ^— 

11 00 
C° = Tt v J ... I dt 1 dt 2 ds 1 ds 2 J dp$°(l3, M 2 /A 2 , m 2 /A 2 ; £ 1; £ 2 , Si, s 2 ) 

00 
xexp [-^°(0,£i,£ 2 ,si,S2)] « 0.3. 

We conclude that the absolute asymptotics of the charge form factor is defined by the 
bubble diagram. The limit Q 2 ^> A 2 of this diagram is due to the ultraviolet regime 
(fl — > 0) in the vertex T M . The vertex is determined directly by the gluon propagator, and 
the 1/Q 2 dependence appears as a manifestation of the ultraviolet behavior of the gluon 
propagator. This is in agreement with the mechanism of hard rescattering and quark 



counting rules [15[]. Namely, the asymptotic behaviour of the form factor is determined 
by the one-gluon exchange between quarks inside a pion. However, in the experimentally 
observed region the triangle diagram dominates in the form factor. 
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V. APPENDIX A 



Below we present a procedure for decomposition of the vertices V a jin written in 
Eq. (p3|). We will consider the vertex for pion with I = n = and = £f = 1/2. 
An extension to the general case is straightforward. As is seen from Eqs. (|10| , [IID we have 
to decompose the operator Fqq in a series over electromagnetic field A. The presence of 
the background gluon field in the vertex is irrelevant to the question under consideration. 
For simplicity we drop the gluon field below. At the end one have to change the partial 
derivatives to the covariant ones. 

The most convenient way is to use the path ordered exponential representation (ei, e 2 
are the charges of u and d quarks): 



oo 



<9 M (x) + z(ei + e 2 )A tl (x) ) /4A^ 
i 



dtPp exp 



t 



4A 2 



d/3 [d (x((3)) + i(ei + e 2 )A(x((3)) 



Pn \ dt 



d A a 



-e o 



1 1 /<-> 

Jdf3a 2 (J3) % J df3a(l3)(d+i(e 1 +e2)A 

n ^ n > 



e o 



7T 



Pr dt da a e 



Using this representation, one gets a nonambiguous decomposition for matrix element of 
the vertex operator: 



e~ ipx F 00 
i 



<9 M (x) +i(e 1 + e 2 )A ll {x)) /4A 



Akx 







dt J da a exp 



—ip x + 



A 



(1 iul 1) | -1 ik \ x I <! >'<'* J ) 
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+i 



r 1 f r 13 

j df3a{(3)A x + ^J df3'a{f3' 
o \ o 



A 



dp'a(p') 



) 



J dt J da a exp | — i(p — k)x — i~j~(p + ^) ^ d(3a{(3) 



x 



1 + i T 6l I d(3a ^ A 



I 



X + 



V 



A 



y d/3'a(/3') 



+ ,^e 2 1 d/3a(/3)A [x - ^ | d/3 , a(/3') J + 0(e 2 ) + O(e^) 

If we now introduce Fourier transformed electromagnetic field, 

d 4 q 



(48) 



(27T) 



then a calculation of coefficients of decomposition to any order in A is reduced to evalu- 
ation of Gaussian path integrals. These integrals can be performed expanding the paths 
a M (/9) into trigonometric series 

oo 

a M (/3) = a° + ^2 [v%cos2irnf3 + v™sin2nnj3} = a° + a M (/3). 

n=0 



In the linear order in A one has to evaluate the integrals: 

1 r /T ~r 1 

(it exp 



rft y <icr a exp 



-i(p — k)x — i—ao(p + k) 



-Hp - k)x - -^(p + kf 



= e- ipx F 00 {-{p + k) 2 )e tkx = e- ipx F 0Q (d 2 (x)) e ikx , 



J da a a\ 



exp 



-i^-a (p + k - (3q) 



= ~i^(p + k- Pq),, exp 



t 

4A 1 



(p + k- (3q) 2 

.\ft .y/i n 

i— a {p + k) + i—paq = exp 



y c?(T a exp ^g— y dP'a(P') 
J da a d^ 



4A 2 



(p + fc - (3qf 



= exp 



4A 2 



P(l-P)q' 



. exp 



- 2/9)g M exp 



4A 2 
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so that 



da a a^((3) exp 
■A 



-i^ip + k) C d(3a((3) + iq^- j ' dffa{& 
A Jo A J 



-*4^(2(p + A:)a»-9a*) ex P 



4A 2 



[(p + A;) 2 -/3(2(p + A;)-g)g] 



Inserting Eq. 



into Eq. (48) one arrives at the representation 

1 2 N 



-ipx 



00 



+(ei - e 2 



9 +i(ei + e 2 )A 

(27T) 



ikx 



e- ipx F 00 (-(p + k) 2 )e ikx 



i i 



^ q+k - p) Mq) J J d pdt—[2{ P + k),-q,} 



o o 



x exp 



t 



■i( P +k) 2 -f3(2(p + k)-q)]q 



-ipx 



4A 2 

Foolf (x)] + (ex - e 2 ) / d 4 yA^y)T^x,y) 



+ 0{e\)+0{e\ 



e ikx + 0(ei) + 0(e 2 2 ) 



with 



r M (s,j/) 



1 1 



(2tt)< 



o o 



dtd(3—[2i d„ (x) -g M ] 



x exp 



4A 2 



9 (x) + /3(2i 9 (x) - q)q 



(49) 



(50) 



To get Eq. ( ^5|) one have to change partial derivatives to covariant derivatives V with the 
vacuum gluon field. 

Equation (|50"D can be written in the form 



e- ipx {F 00 (-(p + kf) + (e 1 -e 2 ) 



X 



(2tt) 4 ~ "^ J 2{p + k)q-q 
F o(-(P + k - qf) - Foo(-(p + A;) 2 )] } e to . 



(51) 



This formula for vertex function coincides with the representations obtained earlier in 
papers ||[|. However, unlike the algorithms used in [|8]|||27| above procedure can be 
applied regularly to calculate the higher order terms, i.e., vertices with two, three or more 
photons. The calculation is reduced to evaluation of Gaussian integrals. 



VI. APPENDIX B 



Here we give some details of calculation of the diagram in Fig 3. The contribution of 
the triangle diagram to the 7177 vertex can be represented in the form (below we use the 
dimensionless, i.e. divided by A, momenta and masses) 
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T% u (p,k h k 2 ) = J da B J - Pl °[ P2 ^ Ps R(h, fc 2 ,p;pi,p 2 ,p 3 ) 



(2vr) ] 



xTrr 3 ij 5 H(p 1 )Q ltl Hip 2 )QjuH(p 3 ). 



(52) 



£ = f d i xd 4 yd 4 ze iipx+klZ - k2y) R{x,y,z;p 1 ,p 2 ,p 3 ) 
R(x,y, z;pi,p2,pa) = exp izBx + ip 3 (z - x 
x exp iyBz + ip 2 (y — 2; 



(53) 



^2 
, V (x) 
F 00 I I oxp 



ix-By + ipi(x — y) 



(54) 



where Q = diag(2/3, — 1/3) is the charge matrix, quark propagator H is defined by 
Eq. (]2?D, and F o is the vertex operator (|rT|). 

Acting by F o onto the exponents in Eq. fl54|) and integrating over the space-time 
coordinates in Eq. fl53"D , we get 



(55) 



lfi7r 4 r r 1 f 

R = ^(p + h + k 2 )—j dtexp y--f3 T M^ - - A {pi+p 3 ) : 



where /^=(/?f/f), 

#i„ = (h +p 2 - pi)ft - ivtf^pt + p 3 )/2, 
/?2/x = (A;i + p 2 - P3V + ivtfpipi + ps)/2, 



Inserting Eq. (|55| ) into Eq. fl5"2|) , one obtains 



?T(P, fei, fe) = ^ (4) (P + h + k 2 )Tr v j^K J da B j dt J 



d 4 p 1 d 4 p 2 d 4 p 3 



(2tt 



1 12 



x exp 



-l^MiP-lipx+pJ 2 



TiT 3 i^H(pi)Q^^H(p 2 ) Qj v H(p z 



(56) 



A calculation of traces gives 



Tf(p, k x , k 2 ) = le^^Xh + k 2 + p) 



xTr v 



3y/2v 



1 1 

IV,. / das J ' ... J dtds\ds 2 ds 3 




1 - sA (\ - s 2 \ (\ - s 3 



1 + SiJ VI + S2/ VI + S3 



1 ml/iv 



d A p 1 d A p 2 d 4 p 3 
(27p 



W a p(pi, P2, P3;m u ) exp [-P T M 2 P + P T Y - h(p, k x ) 
24 



(57) 



where 



W 



a/3 



m u a f3 
^P2P3 



m u a /3 . m u a 

2™ + 1 3PlP2, 

1 Sq 



1 - s 



P T = {pIpIpI) , Y T = (Y?Y?Yf) , 

= t(p + 2A;i) M - if^kj/v, Y 2u = -i(p + 2/c^ + iY^/v, 

= *(p + - if pip - fei)/u, /i = tp 2 /A + ifei + ipfei + ih^f^pv/i 



[2vt + si) / 



^ = 2^ 



-M 2 T 
M 2 r 



-M 2 
(2ut + s 2 ) / 



M 2 
-M 2 



-M 2 T (2vt + s 3 ) I 



M£ v = 2vtS^ - if» v . 



Now one has to integrate over pi, p 2 , P3, and average over directions of the vacuum field 
using the following formulas 

J da B exp(2is 2 4>p~ f f~ f(7 k la ) = J 7 (s 2 if>; kf, Af w ), 

da B e^ a ^f^exp(2ts 2 ^f ya k la ) = -^^(sj^.M,), 
da B ^ vaP ktpnf p exv{2is^p 1 f ia k la ) = 

da B e^pfk{rf K exp(2z S2 # 7 / 7CT A; 1(J ) = e^fcffcf^M; k\, M n ), 



To 



s^ikl + M2) 



sinh 



s 2 v^ 2 + m 2 ; 



cosh 



ip = {siSa + vt(si + s 3 ))/2vx, X = si + s 2 + s 3 + sis 2 s 3 + 2t>t(l + sis 2 + sis 3 + s 2 s 3 ) 



These equations can be derived from the generating formula 



I p2n pit 

— I dip I d9 sin 9 expiiLw Jaw) 

47T JO JO 



sin y 2 (^J^uJfiu i JuuJ) 



au 



2 [JuvJuv i Jui/J/iu 



(5f 



where <£> and 6 1 are the spherical angles, is an anti-symmetric tensor, J Mi , is a dual 
tensor and ± corresponds to the self- and anti-self-dual vacuum field. 

After these manipulations we arrive at the expressions for T£ u (ki, k 2 ,p), Eq. fl28|) , and 
for the transition form factor F^{Q 2 ), Eq. (j36|), with functions $j (i = 1, 2, 3) defined by 
the relations 

$;(M*, Q 2 ; si, s 2 , s 3 , t) = + T x Ri + J^LJ /6t;x 4 - 



The polynomials Pi,Ri, Li have the form 
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Pi = vts\s 2 + (si + vt)[s 1 + s 2 + s 3 + 2nt(l + Sis 3 )], 

P 2 = s 2 [(l + ut(si + s 3 ))(s 2 + 2ut) + si(l + 2wts 3 ) + s 3 (l + 2vtax)\, 

P 3 = vts 2 s 2 3 + (s 3 + i>t)[si + s 2 + s 3 + 2vt(l + S1S3)], 

^ = M*Q u (s, t) - Q 2 Q 2i (s, t) + Qx(s, t), 

Q11 = (si + 2wt)[2sis 3 + wt(si + s 3 - s 2 (l + S1S3))], 

Q12 = 2sis 2 s 3 - 2vt[2 Sl s 3 + s\{l + s 1 s 3 )\ - 2v 2 t 2 {l + sl){s x + s 3 ), 

Ql3 = (>3 + 2vt)[2siS 3 + Vt(s! + S 3 - S 2 (l + S1S3))], 

Q21 = s 2 (si + 2vt)[s 3 + vt(l + sis 3 )} + s 2 (l + 2titsi)[sis 3 + vt(s 1 + s 3 )] 
+vt(si + s 2 )[si + s 3 + 2vt(l + S1S3)], 

Q22 = 2s 2 s 3 + 2wts 2 [si - s 3 + 2s 2 (l + s x s 3 )] + ^^(si + s 3 ), 

Q23 = [«i + s 3 + 2u*(l + sis 3 )][s 2 s 3 + ft(s 2 - s 3 )] + s 2 (s 3 + 2vt)[s 3 + t>t(l + Sis 3 )] 

-s 2 (l + 2wts 3 )[sis 3 + vt(s x + s 3 )], 

Q 31 = 4^x(si + 2vt), 

Q32 = 4v X (s 2 -2vt), 

Q 33 = Av X (s 3 + 2vt), 

Li = s 2 (si + 2wt)[sis 3 + vt(sx + s 3 )], 

£2 = s 2 (s 2 - 2vt)[s 1 s 3 + vtfa + S3)], 

£3 = s 2 (s 3 + 2vt)[s 1 s 3 + vt(s! + S3)]. 

The polynomial x is written in Eq. (|36|). 



As it could be seen above, the background field complicates calculation of diagrams. 
The reason is that in the presence of external field the translation invariance is held only 
for gauge invariant quantities like the diagrams in Fig. 5. This is not true for propagators 
or vertices separately. This means that the energy-momentum is concerved for the total 
diagram but not for the separate vertices and the usual Feynman rules are not valid. So 
that we have to begin calculation in the configuration space. The contribution of the 
triangle diagram to the 7r + 7r~7 vertex has the form 



VII. APPENDIX C 




V (x) 
A 2 



) 



S , (x,y)i7 5 Poo 



/ <->^ 

( v (v) 

A 2 

V 



) 



5(y,^)7M- 



Fourier transformed expression reads 
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h*{k u k 2 ,q) = h*Tr[T+,T-]-Q J da B J 



R 



d A Pid 4 p 2 d 4 p 3 ~ 

^— ^ R(ki,k2,q;Pi,P2,P3) 

d 4 xd 4 yd 4 ze l(klX ~ k2y+qz) R(x, y, z; Pl , p 2 , p 3 ), 



(59) 



R = exp izBx + %p 3 {z — x) 



>2 

V (x) 



"00 



xF, 



A 2 
'v (y) 



00 



V 



A 2 



exp 



exp 



ixBy + ipi(x — y) 
iyBz + ip 2 (y - z) 



(60) 



where Q is the quark charge matrix, H(p) is the quark propagator in the vacuum field 
(P?D, and F 00 is the vertex operator (|TT|). 

The result of action of the operators F 00 and subsequent integration over x, y and z 
looks as 



R = S^{k 1 -k 2 + q) 



where /3 T = 



IQtt 4 



i i 







dtidt 2 
;i + 4t; 2 M 2 ) 2 6XP 



--^FMxP - g(ti,t 2 ;pi,P2,P3) 



(61) 



.Mi 



1 / u(ti + 1 2 )/ Mi \ 

u(l + 4w 2 M 2 ) ^ Ml v(h + At 2 ) I ) ' 
Mf^ = v{tx - 2t 2 )5^ - 2if» v , 

/V = [4(^2 + - Pi)„ + v%t 2 (k 2 + 2pi + 2p 2 + 2p 3 ) fl 
^ivtj^p! + p 3 ) + Uvhf^px + p 2 )]/[4 + u 2 tit 2 ], 

/?2/i = [4(g + p 2 - p 3 )„ + t^Mafy + 2ivt 1 f lt (p 1 + p 3 ) + 2ivt 2 f^ Pl + p 2 )]/[4 + ^ 2 M 2 ] 

# = [*l(Pl +P3) 2 + ^(Pl +P2) 2 + %Vt X t 2 {px + p 2 ) fM f»v(Pl +P3)u]/[4: + v\t 2 ]. 



Using Eq. ( pip in Eq. we arrive at 



1 1 



Aj(fc l5 fc 2 , g) = ^(A* + q - k 2 )Tr v ^ 



dt\dt 2 







do 



B 



d 4 p\d 4 p 2 d 4 p 3 



xTii-f 5 H(p 1 )ij 5 H(p 2 )j^H(p 3 ) exp 



-^p T M!p - g(ti,t 2 ;p 1 ,p2,p 3 ) 



(62) 



As has been discussed above, the delta function in front of the integral in Eq. ( |62|) indicates 
the general conservation of energy-momentum for the total diagram. However, we have 
to evaluate three momentum integrals instead of one loop integral in the usual case of 
zero external field. These Gaussian integrals can be performed as follows. 
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Using representation (27) and calculating the trace of Dirac matrices, one gets 



A£(*i,*2,g) 
l- Sl 



1 - s 2 \ /l - s 3 



X 



1 + s x J U + s 2 



o o 



dt\dt 2 dsids 2 dsz 



l + s 3 

^(px^.ps) exp \-P T M 2 P + P T Y - h(k u q) 



where P T = (pjp^pj), Y T = (y?Y 2 t Y ? 



3 h 



Yf = [v{U + t 2 )k? + v{2t 1 - t 2 )q» - z(l - 2v\t 2 )r(q)]/[v{l + Av\t 2 )], 

Y2 = ~Wi - t 2 )k£ + v(2h + t 2 )<f - 2tv\t 2 r(q) - iffaMvil + Av\t 2 )\, 

Yg = [vih - t 2 )k? + v(2t 1 + t 2 )q» + z(l + 2vH 1 t 2 )r(q) - if" (h)]/[v(l + 4v%t 2 

h(ki, q) = [v(t! + t 2 )k\ + v(Ati + t 2 )q 2 + 2v{2t l - t 2 )qk 1 

+4z(l + vH^qJ^h^/lAvil + Av\t 2 ), 
1 / A x -M 2 M 2 

Mi 



2v(l + AvH.h) 



-Mi A 2 



V Ml 



-No 



-Ni A 3 



A*? = [2v(t! + t 2 ) + si(l + AvHit 2 )\8 i 
= [2v(tt + t 2 ) + s 2 {l + AvH l t 2 )]5 i 
At = [2v{h + t 2 ) + s 3 (l + 4u 2 tit 2 )]5 ( 

= 2u(t! - t 2 )<^ - i(l - Av\t 2 )f^, 
N£ v = 2v(h + t 2 )5^ - z(l + AvH^r, 



and W M has the structure: 



c 



(i-*2)(i-*D 



(1 - s ?)(l- s §) 



^ (l-s?)(l-s|) 



(1 + sis 2 + S1S3 + s 2 s 3 )[^(l - s 2 s 3 ) + if^ u (s 3 - s 2 )], 
(1 + sis 2 + S1S3 + s 2 s 3 )[5^(l + S1S3) + z/^si + s 3 )], 



(1 + S1S2 + S1S3 + s 2 s 3 )[^(l + sis 2 ) + z/m^si + s 2 )], 



£/i,i>aj0 = (<W ~ is 1 f a p)[5f ja ,(s 2 S 3 - 1) + i/^(s 2 - S 3 )] 

+(<f„0 + is 2 f vP )[8 m {l + sis 3 ) + i/^a(si + s 3 )] 
+(5 au + is 3 f au )[S/3^(l + sis 2 ) + ifpn{si + s 2 )]. 
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Integrating over momenta pi, P2 andp3, and, then, averaging over different configurations 
of the vacuum field by means of formulas (k\ = —M 2 , k\ = —M 2 ) 



J da B exp(2#A;i 7 / 7<T g ff ) = Fo^q 2 , M 2 ), 

J da B ifv,{q) exp(2iip fa^f ia q a ) = (fa + k 2 ) ll q 2 T\{^)\ q 2 , M 2 ), 

J daBki a q v qpf^ v f a pexp(2iij)ki^ a q a ) = -(fa + k 2 )^q 2 J 7 2 ('>P; q 2 , M 2 )/2, 



(64) 



Fx 



^q 2 (q 2 + AM 2 ) 
1 



sinh 



cosh 



^q 2 (q 2 + AM 2 ) 



^q 2 (q 2 + AM 2 ) [ 
F 2 = Fo + 4.FiM 

i> = [sis 2 s 3 + vt 1 s 2 (s 1 + s 3 ) + vt 2 s 3 (sx + s 2 ) + vh^sx + s 2 + s 3 + s 1 s 2 s 3 )]/2vx, 
X = 2v(ti + t 2 )(l + sis 2 + sis 3 + s 2 s 3 ) + (1 + 4w 2 tit 2 )(si + s 2 + s 3 + S1S2S3), 



we get Eq. (H) for F*. Functions <$>i(M 2 ,q 2 ; s,t) (i=l,2,3) in Eq. © have the form 
= [F Pi(s, t) + FMM 2 , q 2 ; s, t) + T 2 fa(q 2 - s, t))]/8v X 5 - 



Rj = Q 2 J2 QjkXk, (j = l,3), 

fc=l,2,3 



Q 



Q 



22 



+ 



Q 



2:-! 



R2 ^ [(1 - s\)(l - s 2 ) + (1 - s\)(l - si) ' (1 - - si) 
R 4 = q 2 (q 2 + AM 2 ) £ Q^j + q 2 Qu, 

3=1,2,3 

P J = J2 Qjk^k, (j = 1, 3) 



fc=l,2,3 



21 



+ 



Q 



2r> 



+ 
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-P4 = Q45) Li = L 2 = L 3 = 0, L 4 = — g 2 QikXk- 



fc=l,2,3 



Q11 = (1 + sis 2 + S1S3 + s 2 s 3 )[l + 4w 2 t!t 2 + 2v(ti + t 2 )si + 2v\t 2 (l + 4t; 2 i 1 t 2 )(l - s\)\ 
+ (1 + 4t» 2 t 1 t 2 )(l - sl)[2vtxs 2 (l + sis 3 ) + 2wt 2 s 3 (l + s x s 2 ) + s 2 s 3 - 1] 
+v(t 1 + t 2 )(l-s 2 1 ) X , 

Ql2 = "(I + SlS 2 + S1S3 + S 2 S 3 ) 

2(1 - sl)v 2 t 2 2 (l + AvHj) + (1 + 4v 2 M 2 + 2v(t x + t 2 )s 2 )(l + 2vs 2 t 1 ) 
+ (1 - Sg) 1 + s 2 s 3 + s 3 (st + 8 2 ) + 4utis 2 (l + S1S3) + 11*2(53(1 + sts 2 ) - (si + s 2 )) 
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-Av 2 t x t 2 [s 2 s 3 - 1 + s 3 (s 2 - si) + vt 2 (2s 3 (l + Sis 2 ) + s 2 (l + S1S3) + s x + s 3 )] 
Q21 = (1 + Av\t 2 ) (1 - s?) (1 - si) + v(h - t 2 ) (s 2 (l + Sl s 2 ) + si) 

-(1 + + * 2 )si)(s 2 (si + S 2 ) + 1), 

Q24 = -2x 2 [-^(ti - * 2 )si(s 2 (l + sis 2 ) + si) + + * 2 )(s 2 (si + s 2 ) + 1) 

+Sl(s2(si + S 2 ) + 1 + S1S2)], 

g 23 = (1 + 4u 2 ti* 2 )(l - sl)(l - si) - v(h - t 2 )(s 3 - s 2 )(s 2 s 3 - 1) 

-v{h + * 2 )(s 3 + s 2 )(s 2 s 3 - 1) + 2s 2 (l - s 2 ) + 2s 2 (l - s 2 ) - (1 - s 2 s 3 )(l + s 2 s 3 ), 

g 26 = -2x>(ti - t 2 )(s 2 - s 3 )(s 2 + s 3 ) + v(h + t 2 )(l - s 2 s 3 )(l + s 2 s 3 ) 

+ (s 2 + S 3 )(l - S 2 S 3 )], 

Q 3 i = 2(1 + s x s 2 + S1S3 + s 2 s 3 ) - si) + ut 2 (l + si) + si] 

[Utl(l + Si)+U* 2 (l-Si)+Si], 



g 
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-2(1 - s 2 ) 



>1 + V{U + i 2 ))[s 2 + S3 + + f 2 )(l + S2S3)] 



+V 2 Si(t! -t 2 ) 2 (s 2 + S 3 ) 

g 4 i = (1 + s 1 s 2 + siS3 + s 2 s 3 )[si(l +4u 2 tit 2 ) + + t 2 )(l + s 2 )] 
-(1 + 4v 2 t 1 t 2 )(l - s\) [(1 - s 2 s 3 )(u(ti + t 2 ) + si) + usi(*i - t 2 )(s 3 - s 2 )] , 
g 4 2 = (1 + sis 2 + sis 3 + s 2 s 3 )(s 2 + 2vt 1 )[l + 4v 2 tit 2 + 2v(ti + t 2 )s 2 ] 
+ (1 - s 2 ) s 2 (sxs 3 - 1) + 2vt 1 [s 1 s 2 - 1 + s 3 (s 2 - si)] + 2wt 2 s 3 (si + s 2 ) 

-4w 2 tit 2 [si(s 2 S 3 - 1) + S 3 (siS 2 - 1) + ^l(l + SiS 2 + S1S3 + S 2 S 3 ) 

+vt 2 (s 3 (s 1 + s 2 ) - (1 + SiS 2 ))] , 

g 44 = Avx[d4i + vtid A2 + ^2^43], 

g 45 = -4vx[d 4 4 + ^M 4 5 + ^2^4e], 

d 41 = 4u(*i - t 2 )(s 2 - s 3 )(l + sis 2 + S1S3 + s 2 s 3 ) 

+ (1 - Av 2 tit 2 ) 2(1 + sis 2 + S1S3 + s 2 s 3 )((l + Av 2 tit 2 ){\ + sis 2 + S1S3 + s 2 s 3 ) 

-2(1 - s 2 s 3 )) + (1 + 4t; 2 tit 2 )(l - s\)(l - s 2 )(l - s 2 )" 

+ (1 + 4t; 2 tit 2 )(l - s\){\ - s 2 s 3 )(l + s 2 s 3 ), 

d 42 = 4(1 + sis 2 + S1S3 + s 2 s 3 ) v(ti - t 2 )[si(s 2 - s 3 ) + s 2 s 3 - 1] 

+2(1 - 4u 2 M 2 )s 2 (l + S1S3)] + 2^i(l + Av\t 2 )s 2 (l - s 2 )(l - s 2 3 ), 

d u = As i(l + sis 2 + S1S3 + s 2 s 3 )[v(t 1 - t 2 )(s 3 - s 2 ) + (1 - Av\t 2 )(l - s 2 s 3 )] 

+ (1 + 4t; 2 M 2 )(l - sj)(s 2 + s 3 )(l - s 2 s 3 ) 

d A5 = 4(1 + sis 2 + S1S3 + s 2 s 3 ) v(ti - t 2 )[si(l - s 2 s 3 ) + s 3 - s 2 ] 
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+ (1 - 4^ 2 M 2 )[si(s 3 -s 2 ) + l- s 2 s 3 ]] + (1 + 4v 2 t 1 t 2 )(l - sl)(l + s 2 2 )(l - si), 
Xl = -AvHM 1 - - S3) - s 2 (l - s 2 )(s 2 + 2^0 - s 3 (l - s 2 )(s 3 + 2ut 2 ), 
X 2 = -(1 + 4w 2 tit 2 )(l - s 2 )(l - s 2 ) + 2^(5! + s 3 )(l + sis 3 ) - 2^ 2 s 3 (l - s 2 ) 

+ S 3 (si + S3) + S1S3 + 1, 

V>1 = (1 - s 2 s 3 )[s 2 + S3 + + t 2 )(l + S 2 S 3 )] +f (ti - t 2 )(s 2 - s 3 )(s 2 + S3), 

1p 2 = [1 + S1S3 + S 3 (si + S 3 )](si + + £ 2 )) + - *2)[S1 + S3 + S 3 (l + S1S3)], 

The rest of functions can be obtained by the permutations of arguments t\ <-> t 2 , s 2 <->• S3: 

Ql3(^l? ^2; Si, S 2 , S3) = Ql 2 (t 2 , tl, Si, S3, S2), Q22(^l; ^2; Si, S 2 , S3) = <52l(^2; ^1; Si, S3, s 2 ), 
Q25(^l ? ^2) Si, S2, S3) = Q 2 i(t 2 , tl, Si, S3, S 2 ), Q33 (tl , ^2, Si, S2, S3) = Q$ 2 (t 2 , t\-> s li S3; s 2 ), 
Q 43^1? ^2? Si, S 2 , S3) = Q4 2 (t 2 , ti, Si, S3, s 2 ), 

d^{t X ,t 2 , Si, S 2 , S3) = rf 42 (t2, h, Si, S3, S 2 ), G?46(*l, *2, «1, «2, S3) = d 45 (t 2 , Si, S3, S 2 ), 
X3{h,h, Sl, s 2 , S3) = % 2 (t 2 , ti, Si, S3, S 2 ), -03^1, *2, Si, S 2 , S3) = il> 2 (t 2 , h, Si, S3, s 2 ). 

In the configuration space the 7r + 7r _ 7 vertex corresponding to the bubble diagram in 
Fig. 5 looks as 

A°n(?,y,z) = /i 2 Tr([r + ,r"]_Q) J da B Tr{ij 5 F O0 (x)S(x,y)i-f 5 T fl (z,y)S(y,x)}. 

A calculation of this vertex is analogous to above described case of A^. Using the expres- 
sion for V^z, y) (|25D we get Eq. fl42|) with 

$° = (W, q 2 , M*)P°(s, t) + ^(r; q\ M 2 n )R°(q 2 ; s, t))/{ X °Y- 
P°, R° are the polynomials of t, si, s 2 and S3, 

R°i = <?Y(1 + sis 2 ) 2 (sis 2 + vh( Sl + s 2 ))t 2 /2, 

i? 2 = g 2 [3v(l - 4v 2 tit 2 )(sis 2 + ut x (si + s 2 )) - v 2 tix°]t 2 (l + sis 2 ), 

P 1 ° = X°(l + Sis 2 ) 2 (s 1 -s 2 )t2/2 

P° = 3v{l - Av%t 2 ){s 1 - s 2 )t 2 (l + sis 2 )/2, 



and 



^0 



1 



^y g 2 (g 2 + 4M 2 ) 
1 



sinh 



cosh 



^°vV(9 2 + 4M2) 



^Vg 2 ^ 2 + 4M 2 ) 



^V(? 2 + 4M"2) 
^° = wt 2 /3(sis 2 + u£i(si + s 2 ))/2wx°, 

X° = 2u(ti + t 2 )(l + sis 2 ) + (1 + 4t; 2 tit 2 )(si + s 2 ). 
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TABLES 

TABLE I. Parameters of the model (A 4 = 3B 2 ). 



m u (MeV) 


m d (MeV) 


m s (MeV) 


m c (MeV) 


m b (MeV) 


A (MeV) 


9 


198.3 


198.3 


413 


1650 


4840 


319.5 


9.96 
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TABLE II. The two-photon decay constant g T11 (Gev~ ) and decay 
width r(-7r° — ► 77) (ev); <?* 77 , T* are the values calculated without taking 
into account the spin-field interaction 



(/-7T77 




[30| 


r 


r* 


r cx P | 30 | 


0.235 


0.108 


0.276 


6.3 


1.34 


8.74 
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FIG. 1. A diagram illustrating accuracy of description of meson spectrum. All the masses 
were calculated with the values of parameters given in Table J[ Parameters were fitted using the 
masses of ir, p, K, K* , ifi and T mesons. 
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FIG. 2. A decomposition of propagators and vertices, Eqs. ( p2"tf23| ) into a series in electro- 

(n) 

magnetic field A. A procedure for calculation of the vertices V^- is described in Appendix A. 




FIG. 3. The triangle diagram for the 7*71" — > 7 transition form factor. 



38 



0.30 



0.25 - 



0.20 - 




0.0 2.5 5.0 7.5 10.0 

Q a (GeV 2 ) 

FIG. 4. Transition pion form factor. Solid curve represents a result of the model of induced 
nonlocal quark currents. Long dashed line - calculation without taking into account the spin-field 
interaction. Experimental fit [p^] is given by short dashed line, and Brodsky-Lepage limit EJ 

.186 Gev) is shown by solid straight line; 
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FIG. 5. The triangle (a) and bubble (b) diagrams for pion charge form factor. 
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FIG. 6. The pion charge form factor F 7T (Q 2 ) calculated in the present model (solid line) 
compared with experimental fit (dashed line) 
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FIG. 7. The diagram with /9-meson exchange. 
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